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We consider analytically as well as numerically the finite-size scaling behavior in the stationary 
state near the non-equilibrium phase transition of directed percolation within the mean field regime, 
i.e., above the upper critical dimension. Analogous to equilibrium, usual finite-size scaling is valid 
below the upper critical dimension, whereas it fails above. Performing a momentum analysis of 
associated path integrals we derive modified finite-size scaling forms of the order parameter and 
its higher moments. The results are confirmed by numerical simulations of corresponding high- 
dimensional lattice models. 
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Critical phenomena such as second order phase transi- 
tions are characterized by singularities causing a discon- 
tinuous behavior of various quantities at the transition 
point (e.g. the specific heat, the susceptibility, and the 
correlation length). These singularities are described by 
power-laws defining the well-known critical exponents. 
Studying the phase transition of a given system, the aim 
of investigation is to identify the set of critical exponents 
which characterizes (together with certain universal scal- 
ing functions) the so-called universality class. Since most 
systems are not analytically tractable their critical be- 
havior is often investigated via numerical methods, for 
example Monte Carlo simulations or transfer matrix cal- 
culations. In these cases, the obtained data are limited 
to finite system sizes. Therefore, finite-size scaling (FSS) 
is widely used to extrapolate to the behavior of the in- 
finite systems. In particular, FSS is an efficient method 
to determine the critical exponents and provides certain 
universal scaling functions, i.e., it allows to identify the 
universality class (see 0, H| for reviews) . 

According to the phenomenological FSS theory a 
finite system size L results in a rounding and shifting of 
the singularities. Furthermore, it is assumed that finite- 
size effects are controlled sufficiently close to the critical 
point by the ratio L/£oo- Here, £oo denotes the spatial 
correlation length of the infinite system. Approaching 
the transition point, the correlation length diverges as 
Ceo oc |r — r c | , with the critical exponent v and where 
r is a temperature-like variable which describes the dis- 
tance to the critical point. Finite-size effects are negligi- 
ble for L ^> Coo (i.e. L \t — t c \ v ^> 1). Otherwise, they 
are relevant, i.e., rounding and shifting effects occur if 
L < Coo- It is known in equilibrium that the hypothesis 
of the fundamental role of the ratio L/£oo is valid only 
below the so-called upper critical dimension d c (see |4| 
for a recent review). Above this marginal dimension, 
mean field theories provide exact values of the critical 
exponents as well as of the scaling functions. Rcnormal- 



ization group treatments show that the failure of usual 
FSS within the mean field regime is related to variables 
(scaling fields) which become dangerously irrelevant for 
d > d c 0. Dangerous irrelevant variables affect the 
scaling behavior qualitatively and furthermore cause the 
breakdown of hyperscaling laws which connect the criti- 
cal exponents to the spatial dimension d. Investigations 
of this breakdown of usual finite size scaling trace back 
to 0. After controversial discussions (see e.g. Q and 
reference therein) the problem was recently resolved || 
and a convincing agreement between numerical data and 
field theoretical results is achieved 0, . 

Compared to the equilibrium situation less is known in 
the case of non-equilibrium phase transitions. Therefore, 
we consider in this work the absorbing phase transition 
of directed percolation (DP) as an exemplification. Ac- 
cording to its robustness and ubiquity (including criti- 
cal phenomena in physics, biology, epidemiology, as well 
as catalytic chemical reactions) DP is recognized as the 
paradigm of non-equilibrium phase transitions into ab- 
sorbing states (see [l(j| for a readable review). Although 
an analytical solution is still lacking, DP plays a com- 
parable role in the realm of non-equilibrium phase tran- 
sitions as the Ising model in equilibrium. Previous in- 
vestigations of FSS of DP focus to the absorbing phase 
below d c , where d c — 4 1 lj . Here, we are interested in 
finite-size properties above d c . In particular, we study 
the steady state scaling behavior of finite systems in 
the active phase which is maintained by a homogeneous 
source. Using a momentum space analysis of path in- 
tegrals associated to the field theoretical formulation of 
DP, we derive FSS exponents and universal scaling func- 
tions. Analogous to equilibrium, we demonstrate that 
usual FSS has to be modified in order to describe the 
scaling behavior within the mean field regime. Addi- 
tional numerical simulations confirm the field theoretical 
results. But in contrast to equilibrium we observe a con- 
vincing agreement between the lowest mode finite-size 
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analysis and corresponding numerical results. 

The asymptotic behavior of the DP universality class 
is described by a minimal stochastic Markovian process 
represented by the Langevin equation 0, 0] 



A _1 &n = 



| n- V 2 )n + /i + C (1) 



Here, the density of an active agent n(r,t), defined on 
a mesoscopic (coarse grained) scale, corresponds to the 
order parameter of the non-equilibrium phase transition. 
The control parameter of the transition r attains its crit- 
ical value in an infinite volume at r c . The homogeneous 
source h is conjugated to the order parameter and is usu- 
ally implemented as a spontaneous creation of activity 
(see e.g. 0])- For zero h, a finite positive density occurs 
above the transition point (r < t c ) whereas the absorb- 
ing vacuum state (n = 0) is approached below the transi- 
tion point. Furthermore, £(r, t) denotes the noise which 
accounts for fluctuations of the density n(r, t). This zero- 
mean Gaussian noise represents fast degrees of freedom 
which were eliminated by a suitable coarse graining pro- 
cedure. The noise correlator 



C(r, i)C(r', f ) = A" 1 g> n(r, t) S(r - r') S(t - t') (2) 

is dictated by the existence of the absorbing state n = 0. 

Renormalization group techniques have been applied 
to determine the critical exponents and the universal 
scaling functions [lj El El El El E3 In that case 
path integral formulations are more adequate than the 
Langevin equation approach |l8l llflj ]. Stationary cor- 
relation functions as well as response functions can be 
determined by calculating path integrals with weight 
exp (— J\ where the dynamic functional J describes the 
considered stochastic pro cess. The following dynamic re- 
sponse functional 



jtic proi 
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J = J ^rdtX^h^X^dt + ir-V 2 ) (3) 

+ §(»-»))»-&»} 

is associated to the stochastic process defined by Eq. 
and Eq. J2J. Here, n(r, t) denotes the response field con- 
jugated to the Langevin noise. Furthermore, the coupling 
constants g and g' are equated by an appropriate rescal- 
ing with the redundant parameter K 

n(r, t) -> Jr _1 n(r, t) , n(r, t) -> Kn(r, t) , h -> K h. (4) 

The functional J is invariant under the time inversion 
(so-called rapidity reversal) n(r, t) «-» — n(r, —t) for van- 
ishing (symmetry breaking) source h. 

Using standard techniques known from equilibrium Q , 
it is possible to calculate size-dependent universal scaling 
functions as well as the involved critical exponents. We 
consider DP in a finite cubic geometry of linear size L 



with periodic boundary conditions and expand n and n 
in complex exponential plane waves, e.g. 



n(r,t) =^e^'n(q,t) . 



(5) 



Each component of the wavevector takes only discrete 
values, precisely multiples of 2tt/L including zero. Fol- 
lowing 11], a dynamic free energy functional <£>] for 
the q = mode is constructed by decomposing the crit- 
ical homogenous modes <&(£), §>(t) from their orthogonal 
non-critical complements ^(r, t), \&(r, i), e.g. 



n(r,t) = $(t) + #(r,i) 



(6) 



with $(£) = L d j d d rn(r,t). This leads to a decom- 
position of the response functional J = Jq -\- J\ with 



J(, XL" / di{$(V 1 d t +T+|($-l>))<I>-M>}. (7) 



Now, ^ and \1/ are eliminated by a functional integration 



e -E[*,*] = e -Jb[*,*] 



/x>[#,*]e-^ DM;*.*]. (8) 



The part J\ contributes to the leading scaling behavior 
for d < d c 11]. This will be revisited in a successional 
publication |20(. Here, we consider the mean field regime 
(d > d c ) where J\ provides, besides the shift of the con- 
trol parameter from its mean field value (r" lf = 0) to 
t c x = t c ,oc + 0{L 2 ~ d ) < 0, corrections to the leading 
asymptotic scaling behavior. Hence, we neglect Ji in 
the following but include the shift of the critical point 
to the infinite size value r c = t Ci00 by the redefinition 
t — > t — t c . Correlation functions of the order parame- 
ter (TdUi *(*<*)) =Gfc({ia},r,/i,L d ,A, ff ),canbe de- 
rived from path integrals with weight exp(— £). In that 
way, Eq. J7J and simple dimensional scaling leads to the 
parameter reduction in the correlation function 

G k ({t a }, r, h, L\ A, g) = g- k F k ({Xt a }, r, gh, L d /g 2 ) 
= L- kd ' 2 f k ({L- d ' 2 gXt a }, L d ' 2 T/g, L d h/g) . (9) 

Therefore, it is convenient to define 

<p(s) = L rf/2 $(i) and s^gL^Xt. (10) 

Then Eq. Q yields 

S«Jo= [ ds{<p(d s + T+^(tp-<p)^<p-H(p} (11) 

where the rescaled control parameter and the rescaled 
source 



T = g- x L d / 2 r and H = g^L d h 



(12) 
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are introduced. Note that the whole dynamic functional 
depends on the rescaled parameters T and H only. Fur- 
thermore, the rescaled parameters contain the irrelevant 
parameter g in a dangerous, i.e., singular way. 

The rescaled parameters [Eq. (|12[l ] and Eq. 10 already 
contain a non-trivial result: As usual for critical phenom- 
ena, physical quantities of interest are described in terms 
of generalized homogenous functions. For example, the 
steady state order parameter n — (<£>) and the steady 
state correlation length £ obey for all / > the scaling 
forms (despite of non-universal metric factors) 

n = f' v ' R{Tr 1 / u \hr A / u \Ll), (13) 

e = r x E( T r 1/v *,hr A/l/ ",Li), (u) 

with the universal functions R and S (analogous scaling 



functions are known from equilibrium 4 , 



Usual FSS 



forms involve the correlation length exponent v whereas 
the above modified scaling forms contain the so far un- 
known exponent v* . For d > d c , the order parameter 
exponent j3 and the field exponent A (often called gap 
exponent in equilibrium) equal their mean field values 
[3 = 1 and A = 2. Comparing Eq. |JSJ and Eqs. (|T2|l to 
the scaling forms Eq. and Eq. (|T4*jl for I = L~ x yields 
the FSS exponent for periodic boundary conditions 



2 
d' 



(15) 



Note that v* depends on the spatial dimension in contrast 
to the exponents j3 and A. More important, v* differs 
from the known mean field value of the correlation length 
exponent v = 1/2 for d > d c = 4. Thus, the FSS forms 
are not controlled by the ratio L/£oo °c L\t\ v but by 
L \t\ . This scaling anomaly occurs within the mean 
field regime only and is regarded as breakdown of FSS. 
Furthermore, it can be interpreted as an appearance of 
an additional length scale, termed thermodynamic length 
scale Zoo, which diverges as Zoo <x \ T \~" lilL^J. Similar 
to equilibrium, this length scale coincides with £oo below 
rf c , including v = v* . Thus, the exponent v* fulfills the 
hyperscaling relation v*d = 2/3 + 7' m au dimensions. 

Additionally to the critical exponent v* is it even pos- 
sible to derive universal scaling functions, e.g. i?(0, x, 1). 
The dynamic functional Eq. 11 111 co rresponds to the fol- 
lowing Fokker-Planck equation [ll| 

d s P(tp, s) = [d v [(T + |)v> - H] + 2 |}p(^, a). (16) 

The stationary solution Pq(v?) 



(17) 



can be normalized by an appropriate finite factor C (P, T) 
for H > 0. Straightforward calculations yield the mo- 
ments at bulk criticality (T = 0) 



Thus, the universal FSS functions of the order parameter 
n = ($) = L~ d / 2 (ip), of the order parameter fluctuations 
An = L d (($ 2 ) - ($) 2 ) = (if 2 ) - {ip) 2 as well as of the 



ratios V 

= 1- 



An = 



V 



S = = 1- 



= <<£ 2 >/(<f> 2 

$ 4 )/(3($ 2 ) 2 ) 

2~ r(^) 
L d r(f) 
r^fi) 2 

^r(f) 2 , 



- 1, s = 1 

are given by 



($ 3 )/(3($)($ 2 )), 



(19) 



2r(^±i 

2 ( 1 



2x 




(20) 



(21) 



(22) 
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with the scaling argument x = 2H — 2hL d /g. In contrast 
to equilibrium, the ratios V, S 1 , Q are not finite at the 
critical point (x — * 0). This reflects the different nature 
of the zero order parameter phase in equilibrium and in 
absorbing phase transitions. A ratio that remains finite 
at criticality is obtained via 

2 - 35(a) _ ($ 2 )($ 3 ) - ($)($ 2 )^ _ 1 

~ 2' 



U{x) = 



(23) 



= 2 fc / 2 T(H + k/2) I T(H). 



(18) 



2-3Q(x) ($)($ 4 )-($)($ 2 ) 2 

We expect that this ratio is as useful for absorbing phase 
transition as the Binder cumulant Q is for equilibrium, 
i.e., its value at criticality characterizes the universal- 
ity class. Preliminary numerical investigations below d c 
yield U d=1 = 0.833, U d=2 = 0.704, and U d=3 = 0.61 for 
x -> gg. 

The order parameter n and the ratio Q are shown in 
Fig.ffl Additionally to the above derived universal scal- 
ing functions we plot corresponding simulation results of 
the five-dimensional contact process (CP) as well as of the 
five-dimensional site-directed percolation process (sDP). 
Both models belong to the DP universality class (see [lCJ 
and references therein). In contrast to conventional equi- 
librium simulation techniques, no steady state finite-size 
quantities are available for absorbing phase transitions at 
zero field. Close to the transition point, the systems will 
be soon trapped in the absorbing state without chance of 
escape. As recently pointed out in [23, the natural way 
to circumvent these difficulties is to perform simulations 
in non-zero field at criticality. Thus both, the analytical 
results as well as the numerical simulations reflect that 
well-defined steady state quantities exist close to the crit- 
ical point for h > only. As can be seen in Fig.^ the 
data of the lattice models obey the modified FSS forms 
and the obtained scaling curves are in perfect agreement 
with the results of the continuum theory. A comment is 
worth making: In order to reach numerically the asymp- 
totic scaling regime, the considered system sizes L have 
to exceed all intrinsic non-universal length scales Lq, i.e., 
L Lq. The convincing agreement between the numer- 
ical and the field theoretical results indicates that Lq is 
sufficiently small for the quantities Eqs. (|19I22|) . 
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FIG. 1: The universal order parameter scaling function 
iipbc(0, x, 1) (inset) and the universal fourth order ratio scal- 
ing function Q p bc(0, £, 1) as a function of the rescaled field 
ahh(a,LL) d at criticality for d > d c . The analytically obtained 
scaling functions are in perfect agreement with numerical data 
of the five-dimensional contact process (CP, implemented on 
simple cubic lattices of size L = 4,8,16, A c = 1.13846(11)) 
and of the five- dimensional site-directed percolation process 
(sDP, implemented via the Domany-Kinzel automaton |24| on 
lattices of a generalized bcc-like structure [22J of linear size 
L = 8, 16,32, p c = 0.0359725(2) [||). Note that the numer- 
ical data already belongs to the asymptotic scaling regime. 
In case of the numerical data, non-universal metric factors 
ah and ai, have been introduced in order to norm the uni- 
versal scaling functions, i.e., i? p bc(0, 0, 1) = 1 for the order 
parameter and Q p b c (0, 1, 1) = (bold circle) for the ratio Q 
(see [H|23 for details). 

In contrast to e.g. Q(x), the ratio U exhibits a dif- 
ferent behavior. To be precise, the leading order of U 
is no longer a function of the scaling argument x = 
2hL d /g within the mean field regime. Therefore, non- 
universal corrections to scaling become dominant. Ana- 
lytically, non-universal corrections to U — 1/2 are ob- 
tained by incorporating the shift of the critical value 
(t -► T-OiL^-^l 2 )). The results are confirmed numer- 
ically and will be published in a forthcoming paper |2fj| . 
Again, if the universal leading order of the order param- 
eter moments ($ fc ) is cancelled for ratios such as U , a 
non-universal behavior occurs. Thus U is an appropri- 
ate quantity to investigate the relevance of corrections to 
scaling for d > d c . In summary, a convincing agreement 
is observed between the lowest mode finite-size analysis 
and corresponding numerical results. This is in contrast 
to the situation in equilibrium where it is known that 
the simplest lowest mode approach fails to describe the 
scaling behavior 
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